Translating the COVID-19 epidemiological situation into policies and measures: the Belgian experience.

The COVID-19 pandemic led to sustained surveillance efforts, which made unprecedented volumes and types of data available. In Belgium, these data were used to conduct a targeted and regular assessment of the epidemiological situation. In addition, management tools were developed, incorporating key indicators and thresholds, to define risk levels and offer guidance to policy makers. Categorizing risk into various levels provided a stable framework to monitor the COVID-19 epidemiological situation and allowed for clear communication to authorities. Although translating risk levels into specific public health measures has remained challenging, this experience was foundational for future evaluation of the situation for respiratory infections in general, which, in Belgium, is now based on a management tool combining different data sources.

Outbreak of Corynebacterium diphtheriae among asylum seekers in Belgium in 2022: operational challenges and lessons learnt.

Since 2022, European countries have been facing an outbreak of mainly cutaneous diphtheria caused by toxigenic Corynebacterium diphtheriae among asylum seekers. In Belgium, between 1 March and 31 December 2022, 25 cases of toxigenic C. diphtheriae infection were confirmed among asylum seekers, mostly among young males from Afghanistan. Multi-locus sequence typing showed that most isolates belonged to sequence types 574 or 377, similar to the majority of cases in other European countries. The investigation and management of the outbreak, with many asylum seekers without shelter, required adjustments to case finding, contact tracing and treatment procedures. A test-and-treat centre was organised by non-governmental organisations, the duration of the antimicrobial treatment was shortened to increase compliance, and isolation and contact tracing of cases was not possible. A vaccination centre was opened, and mobile vaccination campaigns were organised to vaccinate a maximum of asylum seekers. No more cases were detected between end December 2022 and May 2023. Unfortunately, though, three cases of respiratory diphtheria, including one death, were reported at the end of June 2023. To prevent future outbreaks, specific attention and sufficient resources should be allocated to this vulnerable population, in Belgium and at international level.

Key risk factors associated with fractal dimension based geographical clustering of COVID-19 data in the Flemish and Brussels region, Belgium.

Introduction: COVID-19 remains a major concern globally. Therefore, it is important to evaluate COVID-19's rapidly changing trends. The fractal dimension has been proposed as a viable method to characterize COVID-19 curves since epidemic data is often subject to considerable heterogeneity. In this study, we aim to investigate the association between various socio-demographic factors and the complexity of the COVID-19 curve as quantified through its fractal dimension. Methods: We collected population indicators data (ethnic composition, socioeconomic status, number of inhabitants, population density, the older adult population proportion, vaccination rate, satisfaction, and trust in the government) at the level of the statistical sector in Belgium. We compared these data with fractal dimension indicators of COVID-19 incidence between 1 January - 31 December 2021 using canonical correlation analysis. Results: Our results showed that these population indicators have a significant association with COVID-19 incidences, with the highest explanatory and predictive power coming from the number of inhabitants, population density, and ethnic composition. Conclusion: It is important to monitor these population indicators during a pandemic, especially when dealing with targeted interventions for a specific population.

Key performance indicators of COVID-19 contact tracing in Belgium from September 2020 to December 2021.

The goal of tracing, testing, and quarantining contacts of infected individuals is to contain the spread of infectious diseases, a strategy widely used during the COVID-19 pandemic. However, limited research exists on the effectiveness of contact tracing, especially with regard to key performance indicators (KPIs), such as the proportion of cases arising from previously identified contacts. In our study, we analyzed contact tracing data from Belgium collected between September 2020 and December 2021 to assess the impact of contact tracing on SARS-CoV-2 transmission and understand its characteristics. Among confirmed cases involved in contact tracing in the Flemish and Brussels-Capital regions, 19.1% were previously identified as close contacts and were aware of prior exposure. These cases, referred to as 'known' to contact tracing operators, reported on average fewer close contacts compared to newly identified individuals (0.80 versus 1.05), resulting in fewer secondary cases (0.23 versus 0.28). Additionally, we calculated the secondary attack rate, representing infections per contact, which was on average lower for the 'known' cases (0.22 versus 0.25) between December 2020 and August 2021. These findings indicate the effectiveness of contact tracing in Belgium in reducing SARS-CoV-2 transmission. Although we were unable to quantify the exact number of prevented cases, our findings emphasize the importance of contact tracing as a public health measure. In addition, contact tracing data provide indications of potential shifts in transmission patterns among different age groups associated with emerging variants of concern and increasing vaccination rates.

Estimated incubation period distributions of mpox using cases from two international European festivals and outbreaks in a club in Berlin, May to June 2022.

BackgroundSince May 2022, an mpox outbreak affecting primarily men who have sex with men (MSM) has occurred in numerous non-endemic countries worldwide. As MSM frequently reported multiple sexual encounters in this outbreak, reliably determining the time of infection is difficult; consequently, estimation of the incubation period is challenging.AimWe aimed to provide valid and precise estimates of the incubation period distribution of mpox by using cases associated with early outbreak settings where infection likely occurred.MethodsColleagues in European countries were invited to provide information on exposure intervals and date of symptom onset for mpox cases who attended a fetish festival in Antwerp, Belgium, a gay pride festival in Gran Canaria, Spain or a particular club in Berlin, Germany, where early mpox outbreaks occurred. Cases of these outbreaks were pooled; doubly censored models using the log-normal, Weibull and Gamma distributions were fitted to estimate the incubation period distribution.ResultsWe included data on 122 laboratory-confirmed cases from 10 European countries. Depending on the distribution used, the median incubation period ranged between 8 and 9 days, with 5th and 95th percentiles ranging from 2 to 3 and from 20 to 23 days, respectively. The shortest interval that included 50% of incubation periods spanned 8 days (4-11 days).ConclusionCurrent public health management of close contacts should consider that in approximately 5% of cases, the incubation period exceeds the commonly used monitoring period of 21 days.

Vaccine effectiveness against transmission of alpha, delta and omicron SARS-COV-2-infection, Belgian contact tracing, 2021-2022.

OBJECTIVES: Vaccine effectiveness against transmission (VET) of SARS-CoV-2-infection can be estimated from secondary attack rates observed during contact tracing. We estimated VET, the vaccine-effect on infectiousness of the index case and susceptibility of the high-risk exposure contact (HREC). METHODS: We fitted RT-PCR-test results from HREC to immunity status (vaccine schedule, prior infection, time since last immunity-conferring event), age, sex, calendar week of sampling, household, background positivity rate and dominant VOC using a multilevel Bayesian regression-model. We included Belgian data collected between January 2021 and January 2022. RESULTS: For primary BNT162b2-vaccination we estimated initial VET at 96% (95%CI 95-97) against Alpha, 87% (95%CI 84-88) against Delta and 31% (95%CI 25-37) against Omicron. Initial VET of booster-vaccination (mRNA primary and booster-vaccination) was 87% (95%CI 86-89) against Delta and 68% (95%CI 65-70) against Omicron. The VET-estimate against Delta and Omicron decreased to 71% (95%CI 64-78) and 55% (95%CI 46-62) respectively, 150-200 days after booster-vaccination. Hybrid immunity, defined as vaccination and documented prior infection, was associated with durable and higher or comparable (by number of antigen exposures) protection against transmission. CONCLUSIONS: While we observed VOC-specific immune-escape, especially by Omicron, and waning over time since immunization, vaccination remained associated with a reduced risk of SARS-CoV-2-transmission.

Fractal dimension based geographical clustering of COVID-19 time series data.

Understanding the local dynamics of COVID-19 transmission calls for an approach that characterizes the incidence curve in a small geographical unit. Given that incidence curves exhibit considerable day-to-day variation, the fractal structure of the time series dynamics is investigated for the Flanders and Brussels Regions of Belgium. For each statistical sector, the smallest administrative geographical entity in Belgium, fractal dimensions of COVID-19 incidence rates, based on rolling time spans of 7, 14, and 21 days were estimated using four different estimators: box-count, Hall-Wood, variogram, and madogram. We found varying patterns of fractal dimensions across time and location. The fractal dimension is further summarized by its mean, variance, and autocorrelation over time. These summary statistics are then used to cluster regions with different incidence rate patterns using k-means clustering. Fractal dimension analysis of COVID-19 incidence thus offers important insight into the past, current, and arguably future evolution of an infectious disease outbreak.

Understanding sexual transmission dynamics and transmission contexts of monkeypox virus: a mixed-methods study of the early outbreak in Belgium (May-June 2022).

OBJECTIVE: The available epidemiological and clinical evidence from the currently ongoing monkeypox (MPX) outbreak in non-endemic areas suggests an important factor of sexual transmission. However, limited information on the behaviour and experiences of individuals with an MPX infection has to date been provided. We aimed to describe the initial phase of the MPX outbreak in Belgium, and to provide a more in-depth description of sexual behaviour and transmission contexts. METHODS: We used routine national surveillance data of 139 confirmed MPX cases with date of symptom onset until 19 June 2022, complemented with 12 semistructured interviews conducted with a subsample of these cases. RESULTS: Sexualised environments, including large festivals and cruising venues for gay men, were the suspected exposure setting for the majority of the cases in the early outbreak phase. In-depth narratives of sexual behaviour support the hypothesis of MPX transmission through close physical contact during sex. Despite awareness of the ongoing MPX outbreak, low self-perceived risk of MPX acquisition and confusing initial signs and symptoms for other STIs or skin conditions delayed early detection of an MPX infection. In addition, we describe relevant contextual factors beyond individual behaviour, related to sexual networks, interpersonal interactions and health systems. Some of these factors may complicate early MPX detection and control efforts. CONCLUSION: Our results highlight the role of sexual contact and networks in the transmission of MPX during the early phase of the outbreak in Belgium. Risk communication messages should consistently and transparently state the predominant sexual transmission potential of MPX virus, and prevention and control measures must be adapted to reflect multilevel factors contributing to MPX transmission risk.

COVID-19 surveillance in the Flemish school system: development of systematic data collection within the public health school system and descriptive analysis of cases reported between October 2020 and June 2021.

BACKGROUND: The age-specific distribution of SARS-CoV-2 cases in schools is not well described. Reported statistics reflect the intensity of community transmission while being shaped by biases from age-dependent testing regimes, as well as effective age-specific interventions. A case surveillance system was introduced within the Flemish school and health-prevention network during the 2020-2021 school year. We present epidemiological data of in-school reported cases in pre-, primary and secondary schools identified by the case surveillance system, in conjunction with test data and community cases from October 2020 to June 2021. METHODS: We describe the development of the surveillance system and provide the number of reported cases and standardized rates per grade over time. We calculated absolute and relative differences in case incidence according to school grade (primary: grades 1-6, and secondary: grades 7-12) using grades 7-8 as a comparator, relating them to non-pharmaceutical infection prevention interventions. Cumulative population incidences (IP) stratified by age, province and socioeconomic status (SES) of the school population are presented with their 95% confidence intervals (CI). RESULTS: A total of 59,996 COVID-19 cases were reported in the school surveillance system, with the highest population adjusted IP in grade 11-12 of 7.39% (95%CI 7.24-7.53) and ranging from 2.23% to 6.25% from pre-school through grade 10. Age-specific reductions in mask introduction and in-person teaching were temporally associated with decreased case incidence, while lower pupil SES was associated with an increase in cumulative cases (excess 2,739/100,000 pupils compared to highest SES tertile). Community testing volumes varied more for children compared to adults, with overall higher child test-positivity. Holidays influence capturing of cases by the system, however efficiency increased to above 75% after further automation and integration in existing structures. CONCLUSION: We demonstrate that effective integration of case surveillance within an electronic school health system is feasible, provides valuable data regarding the evolution of an epidemic among schoolchildren, and is an integral component of public health surveillance and pandemic preparedness. The relationship towards community transmission needs careful evaluation because of age-different testing regimens. In the Flemish region, case incidence within schools exhibited an age gradient that was mitigated through grade-specific interventions, though differences by SES remain.

COVID-19 contact tracing in Belgium: main indicators and performance, January - September 2021.

BACKGROUND: Contact tracing is one of the main public health tools in the control of coronavirus disease 2019 (COVID-19). A centralized contact tracing system was developed in Belgium in 2020. We aim to evaluate the performance and describe the results, between January 01, 2021, and September 30, 2021. The characteristics of COVID-19 cases and the impact of COVID-19 vaccination on testing and tracing are also described. METHODS: We combined laboratory diagnostic test data (molecular and antigen test), vaccination data, and contact tracing data. A descriptive analysis was done to evaluate the performance of contact tracing and describe insights into the epidemiology of COVID-19 by contact tracing. RESULTS: Between January and September 2021, 555.181 COVID-19 cases were reported to the central contact center and 91% were contacted. The average delay between symptom onset and contact tracing initiation was around 5 days, of which 4 days corresponded to pre-testing delay. High-Risk Contacts (HRC) were reported by 49% of the contacted index cases. The mean number of reported HRC was 2.7. In total, 666.869 HRC were reported of which 91% were successfully contacted and 89% of these were tested at least once following the interview. The estimated average secondary attack rate (SAR) among the contacts of the COVID-19 cases who reported at least one contact, was 27% and was significantly higher among household HRC. The proportion of COVID-19 cases who were previously identified as HRC within the central system was 24%. CONCLUSIONS: The contact-tracing system contacted more than 90% of the reported COVID-19 cases and their HRC. This proportion remained stable between January 1 2021 and September 30 2021 despite an increase in cases in March-April 2021. We report high SAR, indicating that through contact tracing a large number of infections were prospectively detected. The system can be further improved by (1) reducing the delay between onset of illness and medical consultation (2) having more exhaustive reporting of HRC by the COVID-19 case.

Vaccine effectiveness against onward transmission of SARS-CoV2-infection by variant of concern and time since vaccination, Belgian contact tracing, 2021.

BACKGROUND: During the first half of 2021, we observed high vaccine effectiveness (VE) against SARS-CoV2-infection. The replacement of the alpha-'variant of concern' (VOC) by the delta-VOC and uncertainty about the time course of immunity called for a re-assessment. METHODS: We estimated VE against transmission of infection (VET) from Belgian contact tracing data for high-risk exposure contacts between 26/01/2021 and 14/12/2021 by susceptibility (VEs) and infectiousness of breakthrough cases (VEi) for a complete schedule of Ad26.COV2.S, ChAdOx1, BNT162b2, mRNA-1273 as well as infection-acquired and hybrid immunity. We used a multilevel Bayesian model and adjusted for personal characteristics (age, sex, household), background exposure, calendar week, VOC and time since immunity conferring-event. FINDINGS: VET-estimates were higher for mRNA-vaccines, over 90%, compared to viral vector vaccines: 66% and 80% for Ad26COV2.S and ChAdOx1 respectively (Alpha, 0-50 days after vaccination). Delta was associated with a 40% increase in odds of transmission and a decrease of VEs (72-64%) and especially of VEi (71-46% for BNT162b2). Infection-acquired and hybrid immunity were less affected by Delta. Waning further reduced VET-estimates: from 81% to 63% for BNT162b2 (Delta, 150-200 days after vaccination). We observed lower initial VEi in the age group 65-84 years (32% vs 46% in the age group 45-64 years for BNT162b2) and faster waning. Hybrid immunity waned slower than vaccine-induced immunity. INTERPRETATION: VEi and VEs-estimates, while remaining significant, were reduced by Delta and waned over time. We observed faster waning in the oldest age group. We should seek to improve vaccine-induced protection in older persons and those vaccinated with viral-vector vaccines.

Gram-negative central line-associated bloodstream infection incidence peak during the summer: a national seasonality cohort study.

Central line-associated bloodstream infections (CLABSI) cause increased morbidity, mortality, and hospital costs that are partially preventable. The phenomenon of seasonality among CLABSI rates has not been fully elucidated, but has implications for accurate surveillance and infection prevention trials. Longitudinal dynamic cohort of hospitals participating in hospital-wide and intensive care unit bloodstream infection surveillance for at least one full year over 2000 to 2014. Mixed-effects negative binomial regression analysis calculated the peak-to-low ratio between months as an adjusted CLABSI incidence rate ratio (IRR) with 95% confidence intervals (CI). Multivariate regression models examined the associations between CLABSI pathogens and ambient temperature and relative humidity. The study population included 104 hospital sites comprising 11,239 CLABSI. Regression analysis identified a hospital-wide increase in total CLABSI during July-August, with a higher gram-negative peak-to-low incidence rate ratio (IRR 2.52 [95% CI 1.92-3.30], p < 0.001) compared to gram-positive bacteria (IRR 1.29 [95% CI 1.11-1.48], p < 0.001). Subgroup analysis replicated this trend for CLABSI diagnosed in the intensive care unit. Only gram-negative CLABSI rates were associated with increased temperature (IRR + 30.3% per 5 °C increase [95% CI 17.3-43.6], p < 0.001) and humidity (IRR + 22.9% per 10% increase [95% CI 7.7-38.3), p < 0.001). The incidence and proportion of gram-negative CLABSI approximately doubled during the summer periods. Ambient temperature and humidity were associated with increases of hospital-acquired gram-negative infections. CLABSI surveillance, preventive intervention trials and epidemiological studies should consider seasonal variation and climatological factors when preparing study designs or interpreting their results.

COVID-19 cases, hospitalizations and deaths in Belgian nursing homes: results of a surveillance conducted between April and December 2020.

BACKGROUND: In Belgium, the first COVID-19 death was reported on 10 March 2020. Nursing home (NH) residents are particularly vulnerable for COVID-19, making it essential to follow-up the spread of COVID-19 in this setting. This manuscript describes the methodology of surveillance and epidemiology of COVID-19 cases, hospitalizations and deaths in Belgian NHs. METHODS: A COVID-19 surveillance in all Belgian NHs (n = 1542) was set up by the regional health authorities and Sciensano. Aggregated data on possible/confirmed COVID-19 cases and hospitalizations and case-based data on deaths were reported by NHs at least once a week. The study period covered April-December 2020. Weekly incidence/prevalence data were calculated per 1000 residents or staff members. RESULTS: This surveillance has been launched within 14 days after the first COVID-19 death in Belgium. Automatic data cleaning was installed using different validation rules. More than 99% of NHs participated at least once, with a median weekly participation rate of 95%. The cumulative incidence of possible/confirmed COVID-19 cases among residents was 206/1000 in the first wave and 367/1000 in the second wave. Most NHs (82%) reported cases in both waves and 74% registered ≥10 possible/confirmed cases among residents at one point in time. In 51% of NHs, at least 10% of staff was absent due to COVID-19 at one point. Between 11 March 2020 and 3 January 2021, 11,329 COVID-19 deaths among NH residents were reported, comprising 57% of all COVID-19 deaths in Belgium in that period. CONCLUSIONS: This surveillance was crucial in mapping COVID-19 in this vulnerable setting and guiding public health interventions, despite limitations of aggregated data and necessary changes in protocol over time. Belgian NHs were severely hit by COVID-19 with many fatal cases. The measure of not allowing visitors, implemented in the beginning of the pandemic, could not avoid the spread of SARS-CoV-2 in the NHs during the first wave. The virus was probably often introduced by staff. Once the virus was introduced, it was difficult to prevent healthcare-associated outbreaks. Although, in contrast to the first wave, personal protective equipment was available in the second wave, again a high number of cases were reported.

Seasonal variation of hospital-acquired bloodstream infections: A national cohort study.

BACKGROUND: Hospital-acquired bloodstream infections (HABSIs) cause increased morbidity, mortality, and hospital costs that are partially preventable. HABSI seasonality has been described for gram-negative bacteria but has not been stratified per infection origin. OBJECTIVE: To assess seasonality among all types of HABSIs and their associations with climate. METHODS: Hospitals performing surveillance for at least 1 full calendar year between 2000 and 2014 were included. Mixed-effects negative binomial regression analysis calculated the peak-to-low monthly ratio as an adjusted HABSI incidence rate ratio (IRR) with 95% confidence intervals (CIs). Another regression model examined associations between HABSI rates and climate variables. These analyses were stratified by microorganism and infectious origin. RESULTS: The study population included 104 hospitals comprising 44,111 HABSIs. Regression analysis identified an incidence rate ratio (IRR) peak in August for gram-negative HABSIs (IRR, 1.59; 95% CI, 1.49-1.71), CLABSIs (IRR, 1.49; 95% CI, 1.30-1.70), and urinary tract HABSI (IRR, 1.52; 95% CI, 1.34-1.74). The gram-negative incidence increased by 13.1% (95% CI, 9.9%-16.4%) for every 5°C increase in temperature. Seasonality was most present among E. coli, K. pneumoniae, E. cloacae, and the nonfermenters. Gram-positive and pulmonary HABSIs did not demonstrate seasonal variation. CONCLUSIONS: Seasonality with summer spikes occurred among gram-negative bacteria, CLABSIs, and urinary tract HABSIs. Higher ambient temperature was associated with gram-negative HABSI rates. The preventable causative factors for seasonality, such as the nurse-to-patient ratio, indoor room temperature or device-utilization, need to be examined to assess areas for improving patient safety.

Establishing an ad hoc COVID-19 mortality surveillance during the first epidemic wave in Belgium, 1 March to 21 June 2020.

BackgroundCOVID-19-related mortality in Belgium has drawn attention for two reasons: its high level, and a good completeness in reporting of deaths. An ad hoc surveillance was established to register COVID-19 death numbers in hospitals, long-term care facilities (LTCF) and the community. Belgium adopted broad inclusion criteria for the COVID-19 death notifications, also including possible cases, resulting in a robust correlation between COVID-19 and all-cause mortality.AimTo document and assess the COVID-19 mortality surveillance in Belgium.MethodsWe described the content and data flows of the registration and we assessed the situation as of 21 June 2020, 103 days after the first death attributable to COVID-19 in Belgium. We calculated the participation rate, the notification delay, the percentage of error detected, and the results of additional investigations.ResultsThe participation rate was 100% for hospitals and 83% for nursing homes. Of all deaths, 85% were recorded within 2 calendar days: 11% within the same day, 41% after 1 day and 33% after 2 days, with a quicker notification in hospitals than in LTCF. Corrections of detected errors reduced the death toll by 5%.ConclusionBelgium implemented a rather complete surveillance of COVID-19 mortality, on account of a rapid investment of the hospitals and LTCF. LTCF could build on past experience of previous surveys and surveillance activities. The adoption of an extended definition of 'COVID-19-related deaths' in a context of limited testing capacity has provided timely information about the severity of the epidemic.

Genomic epidemiology of persistently circulating MDR Shigella sonnei strains associated with men who have sex with men (MSM) in Belgium (2013-19).

OBJECTIVES: Shigella sonnei resistant to first-line antibiotics azithromycin and ciprofloxacin are on the rise globally. The aim of this study was to describe the epidemiology of MDR S. sonnei in Belgium and to identify origins and circulating clusters through WGS. METHODS: We undertook demographic, temporal and geographical analysis of 930 S. sonnei isolates submitted to the Belgian National Reference Centre for Salmonella and Shigella between 2017 and 2019. Phylogenetic analysis of WGS data, genotyping and identification of genetic markers of antimicrobial resistance was performed on 372 Belgian isolates submitted between 2013 and 2019. RESULTS: S. sonnei was identified in 75% (930/1253) of Belgian Shigella isolates submitted between 2017 and 2019. Overall, 7% (69/930) of isolates were resistant to ciprofloxacin alone, 6% (57/930) showed reduced susceptibility to azithromycin alone, and 24% (223/930) exhibited both. Men were at higher risk of carrying a double resistant S. sonnei strain, compared with women (risk ratio = 8.6, 95% CI = 5.4-13.9). Phylogenetic analysis revealed four independent Belgian clusters of persistently circulating MDR strains, associated with men who have sex with men (MSM) and of the same genotypes as previously described international MSM-related clades. Belgian isolates carried various incompatibility (Inc)-type plasmids, the SpA plasmid and ESBL genes. CONCLUSIONS: In Belgium, S. sonnei isolates from men are much more likely to be resistant to important first-line antibiotics than isolates from women. Multiple co-circulating MDR S. sonnei clusters of different genotypes were identified in the MSM community. Further studies on risk groups are needed for targeted prevention, improved clinical and public health management and antimicrobial stewardship in Belgium.

Vaccine effectiveness against infection and onwards transmission of COVID-19: Analysis of Belgian contact tracing data, January-June 2021.

In Belgium, high-risk contacts of an infected person were offered PCR-testing irrespective of their vaccination status. We estimated vaccine effectiveness (VE) against infection and onwards transmission, controlling for previous infections, household-exposure and temporal trends. We included 301,741 tests from 25 January to 24 June 2021. Full-schedule vaccination was associated with significant protection against infection. In addition, mRNA-vaccines reduced onward transmission: VE-estimates increased to >90% when index and contact were fully vaccinated. The small number of viral-vector vaccines included limited interpretability.

Application of a strain-level shotgun metagenomics approach on food samples: resolution of the source of a Salmonella food-borne outbreak.

Food-borne outbreak investigation currently relies on the time-consuming and challenging bacterial isolation from food, to be able to link food-derived strains to more easily obtained isolates from infected people. When no food isolate can be obtained, the source of the outbreak cannot be unambiguously determined. Shotgun metagenomics approaches applied to the food samples could circumvent this need for isolation from the suspected source, but require downstream strain-level data analysis to be able to accurately link to the human isolate. Until now, this approach has not yet been applied outside research settings to analyse real food-borne outbreak samples. In September 2019, a Salmonella outbreak occurred in a hotel school in Bruges, Belgium, affecting over 200 students and teachers. Following standard procedures, the Belgian National Reference Center for human salmonellosis and the National Reference Laboratory for Salmonella in food and feed used conventional analysis based on isolation, serotyping and MLVA (multilocus variable number tandem repeat analysis) comparison, followed by whole-genome sequencing, to confirm the source of the contamination over 2 weeks after receipt of the sample, which was freshly prepared tartar sauce in a meal cooked at the school. Our team used this outbreak as a case study to deliver a proof of concept for a short-read strain-level shotgun metagenomics approach for source tracking. We received two suspect food samples: the full meal and some freshly made tartar sauce served with this meal, requiring the use of raw eggs. After analysis, we could prove, without isolation, that Salmonella was present in both samples, and we obtained an inferred genome of a Salmonella enterica subsp. enterica serovar Enteritidis that could be linked back to the human isolates of the outbreak in a phylogenetic tree. These metagenomics-derived outbreak strains were separated from sporadic cases as well as from another outbreak circulating in Europe at the same time period. This is, to our knowledge, the first Salmonella food-borne outbreak investigation uniquely linking the food source using a metagenomics approach and this in a fast time frame.

Outbreak of Central American born Shigella sonnei in two youth camps in Belgium in the summer of 2019.

In 2019, an outbreak of Shigella sonnei occurred during two youth camps in Belgium. The clustering of isolates from both camps was confirmed by next-generation sequencing, as well as a secondary infection of a technician. The outbreak strain clustered with internationally isolated strains from patients with recent travel history to Central America. This report exemplifies enhanced surveillance and international collaboration between public health institutes by enabling to link local outbreaks to region-specific sublineages circulating abroad.